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Equations are derived which describe the operations (A), for the early stages of the 
run, (J5), after infinite time with no production, and (C), after infinite time at fixed pro- 
duction. The formulas under (A) enable evaluation of the elementary separation 
coefficient, e. Those under (B) enable calculation of the maximum separation in terme of 
various channeling factors which determine the length of theoretical unit cell, h. The 
formulas under (C) enable calculation of optimum production conditions. When measured 
by the smallness of h, the analysis depicts the fractionator as one which is about 100 times 
more powerful than a modern laboratory distillation column. 



I. Introduction 

A mathematical analysis of various factors 
pertinent to the operation of countercurrenl 
electrolysis [1, 2] 3 in capillary tubes, packings, 
and sheets is presented. Special emphasis is 
given to the adverse effects of channeling. Al- 
though the discussion is applied to the specific 
problem of concentrating 39 K, it is applicable to 
open and packed reflux columns in general 
merely by changing the definition of terms. 

The factors involved in the analysis are (1) 
The separation coefficient. This is defined as 
the ratio of the mobilities or velocities of the 39 K 
and 41 K isotopes, e=Vi/v 2 . The actual mechanism 
giving rise to a difference in mobility is not in- 
volved in the discussion. (2) The length of 
column, /?, required to give a separation of e. In 
the discussion the term h takes the place of the 
conventional height of elemental theoretical plate 
or height of theoretical unit. The various ex- 
perimental conditions limiting h are considered 
individually. (3) Factors for continuous produc- 
tion. Calculations based on (1) and (2) are given 



1 This paper will be included in volume 6, division III, of the Manhattan 
Project Technical Series. 

2 Patent Examiner, U. S. Patent Office, and Consultant to Mass Spec- 
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for the most efficient length of tube, the holdup, 
and the yield in grams per kilowatt hour, when 
the product is 39 K/ 41 K=21.3. 
1. Symbols and abbreviations 

A — cross-sectional area of packing, cm 2 
a=lateral spacing between a velocity crest 
and a velocity trough, cm 
Ci,c 2 = concentrations of 39 K, and 41 K, moles/cm 3 
C= normality, moles/liter 
D,D U D 2 = isotropic, lateral and longitudinal co- 
efficients of pure diffusion, cm 2 /sec 
Clover- all or effective longitudinal co- 
efficient of remixing by both pure 
diffusion and convection remixing, 
cm 2 /sec 
€=Vi/v 2 = elementary process separation fac- 
tor 
i^=Faraday cons tant=96, 500 cou- 
lombs/equivalent 
h= length of theoretical unit, as defined in 

text, cm 
rj= viscosity of solution, poise 
I + = total current times positive-ion trans- 
ference number, amperes 
/= length of packing or fractioning path, 

cm 
y= equivalent conductivity of solution 
iVi,iV2=niole fractional concentrations of 39 K 
and 41 K 
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n— dimension] ess exponent 
g=valence of positive ion = l for K 
r,r =radial coordinate and radius of tube, cm 
1?= FJ A Va^ reflux ratio during production 
S t ,S e = over-all separation factor at time t and 
at equilibrium 
tf=time, sec 
v h v 2 = velocity of 39 K and 41 K relative to the 
solution, cm /sec 
Va=mean countervelocity of solution for 
total reflux, cm /sec 
v= local point counter velocity of solution, 
cm /sec 
AF a =mean departure from total reflux mean 
countervelocity, with reference to 
production, cm /sec 
A?;=local point departure from exact ion- 
velocity, solution-velocity balance un- 
der total reflux equilibrium, cm/sec 
V c = cathode volume, cm 3 
x,w= thickness coordinate and semithickness 
between parallel walls, cm 
2=distance from anode face of packing, cm. 

II. Significance of e 

In an electrolytic cell of the type discussed in 
RP1765 [2], the transport mole/cm 2 /sec of the 
light 39 K isotope due to the electric current before 
the liquid counterflow is superimposed is 

CiVi I + 

CiVi+c 2 v 2 qFA' 

in which C\ and c 2 are mole concentrations of 39 K 
and 41 K per milliliter at time t; v L and v 2 are their 
velocities; I + is the positive-ion current equal to 
the total current times the cation transference 
number; 2^=96, 500 is the Faraday constant; A is 
the effective cross-sectional area; and q is the 
valence of the positive ion, being unity for K+. 
When liquid counterflow is superimposed at a 
rate intermediate between the velocities of the 
39 K and 41 K, these are transported in opposite 
directions and an isotope concentration gradient 
builds up. This gradient begins adjacent to the 
cathode and works its way to the anode, finally 
reaching an equilibrium value. When this equi- 
librium is reached, the stream flow plus the back 
diffusion and convection remixing just balances 
the 39 K ion current. If the potential is then re- 
moved, the flow of 39 K in the anode direction is 



V a c 1 -\-D'(dci/dz) t = a > moles per square centi- 
meters a second, where V a is the mean liquid 
velocity, D f is the over-all or effective coefficient 
of remixing by both diffusion and convection 
processes, and z refers to length along the tube. 
When the v average K ion is static against the 
moving stream, as is the case during perfect 
operation, and when equilibrium is reached for 
the entire tube, conditions in the system are 
expressed by 



Cltfi 



/+ 



Ci#i-fc 2 ^2 qJ?A 



^Vad + D'idcjdz)^ 



(1) 



where (dci/dz) t=co represents the maximum at- 
tainable gradient. 

The liquid countervelocity, V aj is so adjusted at 
the cathode end that the net forward transport of 
positive ions by the current is balanced by the 
backward transport of the stream, so that 
I + l(qFA)=--V a (ci+c 2 ). Hence there is no net 
positive-ion current after this adjustment. The 
small currents of 39 K and 41 K that cause the tube 
to build up a gradient are equal and opposite, and 
finally disappear when the equilibrium gradient is 
reached. Hence, under proper conditions, the 
entire current is carried by the negative ions 
throughout the run. V a is maintained constant 
throughout the run and has the value 



v a - 



~qFA( Cl + C 2 ) 



Substituting eq 2 in oq 1 and collecting 

/ + C\C 2 (Vi — V 2 ) 



qFA (ci+c 2 )(c l Vi + c 2 v 2 y 



--D' '(dcjdz) ,. 



(2) 



(3) 



As V\ and v 2 are nearly equal, it is sufficiently 
accurate to write (ciVi J rC 2 v 2 ) = (ci J rc 2 )v 2 and, also, 
to express isotope concentrations as mole fractions 

by Ni = Ci/(ci + c 2 ), N 2 = c 2 l{c l J rC 2 ), and dc x = 
(ci + c 2 )dAfi- With these relations, eq 3 becomes 



I+^-v,) 



dz~- 



dN x 



qFAD'( Cl +c 2 )v 2 N X N< 



(4) 



As JVi + #2= 1, and hence dN x = —dN 2) it follows 
that 

dN 1 /N 1 N 2 =dN 1 /N l -dN 2 /N 2 =d ln(iVyiV 2 ), 

and the integration of eq 4 over the length, I, of the 
tube, assuming uniform salt concentration, gives 
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I+ivi-v^l 



In 



M/tf,).-, 



2 ^£»'( Cl +c 2 )t> 2 (iV,/^) 



(5) 



This is the equilibrium equation for the tube and 
applies only after the gradient ceases to rise in 
any part of the tube, 

The elementary process separation coefficient 
is defined as the ratio of the mobilities or velocities, 
of the light to the heavy-isotope ions, thus 



e=v 1 fv 2 or €— 1 = 



V1 — V2 



v 2 

The over-all equilibrium separation factor is 
defined as the ratio of the abundance ratio at the 
cathode to the abundance ratio at the anode, thus 






(6) 



As * is close to unity, e— 1 can be replaced by 
In 6. Equation 5 therefore can be written 



In S.< 



I+(e-l)l 



I+l In e 



or 



iFAD'id+Ci^qFAD'id+Ci) 

S e =(e)"", 



(7) 
(8) 



where h is the height of theoretical unit and equals 
qFAD'(c l -\-c 2 )II+ = D'/Y in a ad l/h is the number of 
theoretical units contained in the length I. 

All quantities are measurable in eq 7 and 8 
with the exception of the over-all remixing' coeffi- 
cient D'; this can be calculated under certain 
simplified experimental conditions, for example, 
those prevailing in a single vertical capillary tube. 
Equations 7 and 8 thus make it possible to calcu- 
late the maximum separation attainable. 

At the beginning of the run there is no remixing 
as dci/dz=0, and the total transport of light iso- 
tope is concentrated in the cathode compartment. 
The initial transport, in moles of 39 K per second, 
from eq 3 is 



7+ N.N^-v,) 
qF V 2 



7+ Ayy 2 (<-i) 

qF 



(9) 



The initial rate of rise in concentration at the 
cathode is inversely proportional to the cathode 
volume. When ANi is the mole fraction gain in 
concentration of the light isotope in t seconds; 
V c , the cathode volume; and C, the normal- 
ity; the moles of light isotope accumulated is 
VcCAN^Q- 3 moles 39 K. This is equal to the rate 
given by eq 9 times the time. Hence 



VcCANxlO- 



qF 



or 



VJJqF AN, 
lOOOl+tNA 7 ? 



(10) 



As AjVx is small compared to \ { during initial 
operation, 

£^Aln(NJNi) =hi(N 1 /N i ) t - 

ln(M/Ay ( =o=lnS<. 

Equation 10 now becomes 



Vj OgFinS, 

IlKlii/ / 



(11) 



With eq 11 it is possible to calculate the ratio of 
mobilities from the initial slope of the curve of 
In S t plott ed wit h respect to t. If S t is near unity, 
it is sufficiently accurate to write In S t = S t —l 
and plot this w i t h respect to t. Also Ni/N 2 may be 
plotted, but in this case the initial slope must be 
divided by (Ni/N 2 ) (mean) in making the calcu- 
lation of e. 

Equation 10 or 1 1 applies regardless of whether 
the light isotope is rare or abundant. In cases in- 
volving the concentration of the heavy isotope, 
.Y, and jV 2 are interchanged, and S t is defined as 
(iVi/iVi) ^/(JVi/iVi) f -o- The definition of e then 
remains the same, namely, the ratio of the light 
to the heavy-ion mobilities. 

III. Calculation of h 

The factor h is defined for the present process 
as that segmental length of tube or cell that will 
boost the abundance ratio by a factor e when the 
concentration gradient has built up to a steady 
state under " total reflux", i.e., cathode sampling 
rate is zero. It is important that h be kept small 
as the power requirement, holdup, and time to 
reach equilibrium are each proportional to it. 
This section discusses conditions under which a 
small h can be realized, as well as other conditions 
that will lead to a large h. 

According to eq 2, 7, and 8, the length h is 
given by 

h=D'/V a) (12) 

where D' is the over-all longitudinal remixing 
coefficient, and V a is the mean liquid counter- 
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velocity that balances the mean drift velocity 
of the K ions. In all cases met in practice the 
effective coefficient D' is made up of both pure 
molecular back diffusion and convective remixing. 
The latter is due to convection of K resulting from 
the fact that the point liquid countervelocity is 
not in exact balance with the point ion velocity 
over the entire flow cross section, being higher 
than the ion velocity in one region and lower than 
the ion velocity in another. It is necessary to 
evaluate this latter remixing effect to arrive at a 
value for D' or h. As convective remixing depends 
not only on the longitudinal velocity of the K 
but is diminished by the cross diffusion between 
the counterflowing streams of K, it is necessary 
to analyze the steady state in terms of the profiles 
of longitudinal velocity and lateral diffusivity. 
This can be done rigorously only for cases involv- 
ing geometrical symmetry, in a manner now to be 
discussed. 

Consider a round tube that may be either a 
single narrow capillary or a much larger packed 
tube. Assume that the mean liquid counter- 
velocity is held in exact balance with the mean ion 
velocity, notwithstanding that the center of the 
tube has a liquid velocity greater, and the periph- 
ery a liquid velocity less, than the ion velocity 
in these regions. Assume further that the profile 
of the velocity departure is symmetrical to the 
center of the tube. Then, for the point difference 
between liquid countervelocity, v, and 39 K ion 
velocity v x 

kv=Vi—v, (13) 

where the point velocity departure, Av, is a sym- 
metrical function of r. 

Also, as there is a net balance of transport with 
respect to average values, 



v l c 1 +v 2 c 2 = V a (c 1 +c 2 ), 



(14) 



where v { and v 2 are the averaged point 41 K ion 
velocities and the c's are the averaged point 
concentrations of the two isotopes. 

As by definition, e=Vi/v 2 , which is closely 
Vi/v 2 , and as e— 1 is quite small compared to unity, 
eq 14 may be written 



Th~V a =V a c 2 (e-l)/(c 1 + c 2 ). 



(15) 



integral equations that describe the steady state 
operation of the tube. 

Let dcjdz be the longitudinal concentration 
gradient finally established after very prolonged 
running, c x being the average 39 K concentration 
over the cross section at level z. This gradient is 
also the gradient for each point in the cross section 
because the distribution curves for concentration 
along the wall, along the axis, or along an inter- 
mediate cylinder, are closely parallel in all cases 
to be considered. 

Consider an element of volume at a point r, z. 
Due to such lack of equality as may exist between 
ion and liquid velocity, the 39 K is removed longi- 
tudinal from the element at a rate {i\ — v)dci/dz 
moles sec -1 cm" 3 . But as point concentrations 
in the tube have reached a steady state, this 
removal must be exactly compensated by a lateral 
diffusion of 39 K into the element at a rate (d/rdr) 
(D x r dcjdr) moles sec" 1 cm" 3 , in which D x is the 
radial coefficient of diffusion. Hence the steady 
state requires 

djD^ddfdr) = (v x - v) dcjdz. (1 6) 

rdr 

A second equilibrium condition, assuming zero 
sampling rate at the cathode, is that the integrated 
transport of 39 K is zero. This is expressed by 



/; 



\(vi-v)c 1 -D 2 dc 1 ldz]2wrdr=0 y (17) 



Equations 13 and 15 will be returned to, it now 
being necessary to consider the differential and 



where (i\ — i^)Ci is the moles/cm 2 /sec transported at 
the point r,z due to the velocity departure, and 
D 2 i>Ci/bz is the moles cm" 2 sec" 1 diffusing back, 
D 2 being the longitudinal coefficient of pure dif- 
fusion. 

Equations 16 and 17 can be combined after 
suitable transformation as follows: Integration of 
eq 16 gives 

dc l = c)c l /dz{~ f 'fa— v)rdr\drjD x r } (18) 

whereas integration of eq 17 by parts gives 

Ci I °(v 1 — v)2Trdr— ° i (v l —v)27rrdrdc 1 = 
Jo Jo Jo 

wrlD.dcJdz. (19) 



For dci in eq 19 is substituted the expression given 
by eq 18, giving 
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Ci r°(v 1 -v)2rdr=l2 f °[ C^-nn/rl 

drjDxr+DA&Jbz. (20) 
But 

°v i 2rdr=v 1 ro and I °v2rdr=V a rl. 

Jo Jo 

Hence returning to eq 13 and 15, and combining 
these with eq 20, 



=I 2 I'°[J>H 

drjD l r+D 2 rl\(U l ld^ (21) 



'c&Vaje— \)r 2 Q 

Ci + C 2 



Now let dz=h. Then by definition of h the 
corresponding longitudinal increase in the abund- 
ance ratio is by a factor e, where 

c 2 +dc 2 

As dci=—dc 2 , and is small, this last expression 
can be written 

(gi+~2)<fci) 

c x c 2 



I 



and substitution in eq 21 gives 
V a rl 



+D 2 /V a . (22) 



Equation 22 is the general result sought. It 
enables the calculation of h for any round tube 
whether packed or open, provided Av and D t are 
known symmetrical functions of the radial coor- 
dinate, r. 

In a similar fashion it can be shown that h for 
a long, thin packing or film, contained between 
parallel walls, where flow is lengthwise and a 
uniform balance can be assumed breadthwise, is 
given by 



r wr rx ~12 

I (Av)dx dx/Di 
h= JoUo ^, - + D 2 /V a) 



VaW 



(22a) 



where Av is the symmetrical velocity departure as 
a function of the thickness coordinate, x, the 
origin of the latter being midway of the walls, and 
w is the semispacing. In general, the upper and 
lower limits of x should be chosen such that 
dci/dx=0 at both limits. 

Five special cases under eq 22 and two under 
eq 22a will now be considered in detail. 



Case I. The Isothermal Capillary. — Consider a 

single narrow capillary tube and assume an inap- 
preciable rise in axial temperature relative to the 
wall tempera tun?. Then the ion velocity is uniform 
over the cross section so that v x = Vi. But the 
liquid countervelocity is a symmetrical parabolic 
function of r, and this is true even if the flow is 
aided or impeded by electro-osmosis, as it is 
known that the latter type of flow also has a 
parabolic profile [3]. 
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Figure 1.- — Profile of isothermal capillary tube. 

The profile of the velocity departure is illus- 
trated in figure 1 and is given by 

Av=Vl - v==zVl - v= V a [2(r/r ) 2 -l]. (23) 

This is to be substituted in eq 22. The profile of 
the radial diffusion coefficient D x is a horizontal 
line in this case, and hence A may be brought 
outside the integral sign. The diffusion coeffi- 
cient vertically will be the same as radially, and 
both will be designated I). 

For the first operation indicated by eq 22, 
f J Q f (At;)rdr] = (F a 2 rJ/4)[(r/r ) 8 -2(r/r ) 6 + (r/r ) 4 ]- 
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The remaining integration with limits r , and 
leads to 

h=(V a ri/2D)(l/8-l/Z + l/4)+DIV a 
= (l/48)V a Ti/D+DIV a . (24) 

From this result it is seen that the minimum value 
of h is r /2-y3, and that the liquid countervelocity 
to achieve this minimum is 

V a =(D/r )4^. 

The values of h for a 1.0-mm inside diameter 
isothermal capillary for several voltage gradients 
lengthwise of the tube are listed in table 1. The 
temperature is assumed to be 18° C throughout, 
for which the mobility in 1.0 JV KC1 [6] is 
0.000498 cm sec" 1 volt" 1 . The value of D in 
KC1 solution up to 1.0 JV is close to 1.55 X10" 5 
cm 2 /sec. 

Table 1. — Height of theoretical unit of 1.0-mm inside 
diameter isothermal capillary at 18° C 



Potential 


Va 


(1/48) V a rllD 


D/Va 


h 


ijcm 
1__ 


cm/sec 

0. 000468 

. 00215 

. C04C8 


0. 00168 
. 00722 
.0168 


0.03110 
. 00722 
.00311 


cm 
0. 0328 
.0144 min. 
.0199 


4.31... 


10-. 





In practice, because of PR heating, the assumed 
isothermal condition is not realized. If the capil- 
lary is set vertically with cathode feed at the top, 
and say a 40° C axial rise in temperature is per- 
mitted, it is possible to so select the inside diameter 
that thermal convection up the center and down 
along the walls will partially balance the excess 
and deficiency of the countervelocity. The re- 
sultant profile, illustrated in figure 2, is a closer 
approach to ideal point balance than figure 1. 
In consequence, h will be less than the values 
indicated in table 1. Calculations indicate that 
a vertical capillary under perfect control with an 
inside diameter of 0.7 mm, and a gradient of 15 
v/cm, can yield an h of 0.002 cm, or 500 theoretical 
units per centimeter length. 

If the capillary is set horizontally, the effect of 
heating results in a lateral convective mixing. 
This increases the value of D in the first term of 
eq 24. Hence the value of h will be somewhat 
lower than indicated in table 1, provided the in- 
creased back diffusion represented by the term 
D/V a is unimportant, as is the case for voltage 
gradients higher than about 10. 



The foregoing calculations apply only to a single 
capillary tube. To be applicable to a bundle of 
tubes it would be necessary for the mean ion and 
liquid velocities to balance in each. The exact- 
ness of this individual tube balance will be realized 
when it is considered that maximum separation 
can occur in each tube only when the mean liquid 
velocity has a value lying between the velocities 
of the 39 K and 41 K ions. As, for a given liquid 
head, the liquid velocity varies as the square of 
capillary bore, and as the two ions differ in veloc- 
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Figure 2. — Profile improved by thermal convection. 

ity by only about 0.3 percent, a suitable bundle 
must consist of tubes with variations in bore of 
less than 0.2 percent. Such a requirement can- 
not be met by selection of ordinary glass tubing, 
although extruded tubes might have the required 
uniformity. However, a further requirement is 
that the tubes be at nearly the same temperature, 
and this is difficult to achieve as the solution has 
a negative temperature coefficient of resistance. 

In contrast, if the above bundle is provided 
with numerous lateral communications between 
the tubes, the resulting cross diffusion will greatly 
compensate for lack of individual tube balance, 
and it becomes possible to achieve a reasonably 
low h. This action, though on a larger lateral 
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scale, is similar to the effect of radial diffusion in 
the single capillary, which partially compensates 
for the point lack of balance over the cross section 
and gave rise to the reasonably small first term 
in eq 24. 

A tube packed with sand, beads, etc., is analo- 
gous to a capillary tube bundle having cross 
communication between passages. It w r ill be 
clear, therefore, that isotope separation in other 
than an ideal packing is possible only because of 
cross diffusion. Various cases of nonideal pack- 
ings will be treated hereinafter, but first it appears 
desirable to discuss the ideal, though as yet non- 
existent, packing in order to show the ultimate 
possibilities of this method of fractionation. 

Case II. The Ideal Packing. — The perfect pack- 
ing should have the properties of an ideal tube 
bundle. This requires that the pores be small 
and nearly identical in shape, cross-sectional size, 
and length. The ideal case must exclude the 
presence of a temperature gradient normal to the 
direction of liquid counterflow because fluidity of 
the solution increases with a somewhat greater 
temperature coefficient than does the mobility of 
the ions. This last consideration suggests that 
the ideal packing must be a face-cooled diaphragm. 

In order to maintain maximum temperatures 
below 100° C, calculations indicate that the dia- 
phragm should have a thickness not exceeding 
about 1 cm. Where one diaphragm is insufficient 
to give the desired separation factor, a series can 
be used with intermediate spaces for insertion of 
cooling elements. 

A further requirement is that thermal convection 
across the thickness of the diaphragm must be 
eliminated. Consider a frit made from 50-mesh 
beads. With the faces vertical, a convection 
current will occur inwardly from both faces at 
the bottom, up the midsection, and outwardly 
toward both faces at the top. The convection 
at top and bottom will cause remixing. On the 
other hand, with the diaphragm horizontal, the 
upper half will comprise cool liquid overlying hot 
and the thermal instability will result in convection 
currents toward and away from the upper lace. 
Calculations indicate that if thermal convections 
are to be reduced to negligible velocity the frit 
should be made from beads of about 200 mesh 
per inch or finer. This corresponds to a mean 
pore radius of 0.002 cm or less. 

In summary, it appears that the ideal packing 



should be a face cooled diaphragm of extremely 
uniform thickness and porosity with mean pore 
radius of less than 0.002 cm. For such a dia- 
phragm the first term of eq 24, using for r the 
mean pore radius, is negligible; hence, h=D/V ay 
where D is the longitudinal diffusion coefficient 
of the free liquid and V a the mean liquid counter- 
velocity in the pore. The following are the 
values of h under the ideal conditions assumed 
when V a is taken as 0.00498 cm/sec, corresponding 
to a 1.0 N KC1 solution and an applied gradient 
of 10 v/cm at 18° C: 

Values of D at temperatures other than 18° C 
were estimated by assuming proportionality to 
absolute temperature times the fluidity of the 
water for that temperature [81. 

At 18° C; /}=1.55X10~ 5 cm 2 /sec; A=0.0031cm. 

At40°C; Z?=2.68X10- 5 cm 2 /sec; /^ = 0.0054cm. 

At80°C; 77=5.57X10- 5 cm 2 /sec; h=().0U2cm. 

The h at 18° C would apply to diaphragm 
regions close to either face. The // at 80° C w T ould 
be encountered across the midplane of a frit 
having this maximum temperature. As the 
mobility increases toward the midplane, the volt- 
age gradient across this plane will drop to a value 
such that the ion velocity and liquid counter- 
velocity are in balance. 

Case III. Packed Tube With Unijorm Poros- 
ity and Hot Core. — Assume first that the packing 
has a pore radius of about 0.002 cm or less, so 
that thermal convection is negligible. The value 
of h will be DjV a at the cooled end faces of the 
packing just as in case IL At points remote 
from the ends, however, the ion velocity and liquid 
velocity are not in point balance over the cross 
section due to the fact that the temperature coef- 
ficient for fluidity is somewhat higher than for 
mobility. The bulk profile of the velocity de- 
parture will approximate the appearance shown 
in figure 3, and will be assumed to be parabolic, 
as given by the relation 



Av=Av m [2(r/r ) 2 -H 



(25) 



where all velocities and diffusion coefficients are 
now to be considered as superficial. The super- 
ficial values are obtained by multiplying the mean 
pore values by the fraction voids. 
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Figure 3. — Profile of 'packed tube, uniform porosity, hot 
core. 

Equation 22 is applicable. A rigorous solution 
requires that the lateral diffusivity D x be expressed 
in terms of r, as D x increases with temperature 
and has its maximum value at the point of maxi- 
mum temperature, r=0. It is sufficiently ac- 
curate for our purpose here to consider D x as 
having a uniform value based on the mean tem- 
perature between wall and axis. The result of 
the integration is then 

A=(l/48) (AvJVayVarl/Dt+Dz/Va, (26) 

where V a is equal to the superficial mean ion 
velocity based on the mean temperature. 

As eq 26 involves ratios of velocity to diffusivity, 
each of which is to multiplied by the fraction voids 
to obtain superficial values, it is just as proper to 
use mean pore values. Further, as V a is adjusted 
to balance the mean ion velocity, and as the 
latter is closely the same function of temperature 
as is the diffusivity, calculations can be based on 
KC1 solution at 18°, for which, at 1.0 normal, 
F a =0.000498 cm sec" 1 volt" 1 , and D^D 2 = 
1.55X10- 5 cm 2 /sec. 

Consider a tube packed with 200-mesh sand and 
having an applied gradient of 10 v/cm. The 
maximum permissible temperature rise is assumed 
to be 80° C. For this rise, due to the fact that 
fluidity increases more rapidly with temperature 
than mobility, the ratio Av m /V a is calculated to be 



0.10. Table 2 lists the values of h for various 
tube sizes. 

Tablf 2. — Packed tube with uniform porosity, 80° C 
temperature (rise, 10 v/cm) 



Inside diame- 
ter of tube 


h 


cm 


cm 


0.5 


0. 0073 


1.0 


.0198 


2.0 


.0701 


4.0 


.271 


8.0 


1.07 



It should be emphasized here that the values in 
table 2 are for an extreme case, remote from the 
ends of the packing, and with no thermal convec- 
tion. If the tube is set vertically, and instead of 
200-mesh sand, a coarser packing of about 100 
mesh is used, calculations indicate that thermal 
convection will largely overcome the lack of point 
balance indicated in figure 3. If a tube of 100- 
mesh sand is set horizontally, thermal convection 
will be sufficient to markedly increase the effective 
lateral diffusion coefficient, D\. For these reasons 
it is possible bo obtain lower values for h than are 
indicated by table 2. On the other hand, if the 
packing contains gas bubbles, the effective value 
of Di is decreased; hence, h will be increased. 

Case IV. Packed Tube With Wall Channeling. — 
Uniform porosity was assumed in case III. 
Actually, when a tube is packed with beads or 
sand, the porosity adjacent the wall, and for a 
width of about half a bead, is about 50 percent 
greater than elsewhere. This fact can usually be 
disregarded when the beads are very small and 
held under pressure. In the case of frits the 
higher porosity at the wall can be eliminated 
with suitable cement. It is of interest, however, 
to calculate the value of h for cases in which, due 
to failure of one kind or another, a wall channel 
makes its appearance. 

Consider first a wall channel of uniform thick- 
ness over the entire circumference. A suitable 
function for Av is 



A#= A?; m — (Av w +Av m ) (r/r Q ) 2 



(27) 



where Av w is the maximum velocity departure 
near the wall. The parameter, n, determines the 
extent of the flat portion of the velocity profile 
and for large values of this parameter the profile 
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is flat over the major portion of the cross section, 
as illustrated in figure 4. 
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Figure 4. — Packed tube with symmetrical wall channeling. 

Equation 27 must satisfy the condition for 
average balance, 

fro 

(Av)2Trdr=0, 

from which it is found thai 

n=Av w /Av m . 

From eq 22, neglecting the D 2 /V a term, it is 
found that 

h=(Avy /2V a D 1 )[l/4-l/(n+2) + l/Hn+l)] f (28) 

and for those cases in which the wall velocity 
departure is 10 or more times the departure in 
the packing, this becomes 



h^0.1(Av m rolV a ) 2 V a ID lt 



(29) 



In terms of the channel thickness, d, this can 
be written 



It is seen from eq 29 that the values of h will be 
about five times those for case III for the same 
values of Av m , thus showing the adverse effect of 
abrupt as compared with gradual departure from 
point balance. A symmetrical wall channeling 
sufficient to cause a velocity departure in the 
packing of 10 percent will lead to a value of h of 
0.12 cm for a 1.0-cm inside diameter tube oper- 
ating with an applied voltage gradient of 10 v/cm. 

From eq 29a it is seen that if the thickness of 
the wall channel is the same, the h for various 
tube sizes will remain the same. This is because 
the wall-velocity departure for a given head de- 
pends only on the channel thickness. 
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Figure 5. — Packed tube — wn symmetrical wall channeling. 

Although a rigorous treatment cannot be given 
for cases involving unsymmetrical channeling, 
such as illustrated in figure 5, it can be stated 
that for the same Av m the value of h will be much 
higher than in the symmetrical case. Thus, in 
the general case, h varies as the square of the 
average diffusion distance between the channel 
and the interior of the packing and inversely as 
the lateral diffusion area between the channel and 
the packing. Hence, a wall channel due to a 
localized breakthrough may easily result in such 
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a large h that isotope separation for a 10-cm 
length of packing will be inappreciable. 

Case V. Core Channeling. — Assume the forma- 
tion of a single narrow round channel of radius 
r x and coaxial with the packed tube. The profile 
of the velocity departure will be parabolic in the 
channel and essentially uniform in the surround- 
ing packing, as illustrated in figure 6. Hence 
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Figure 6. — Packed tube — symmetrical core channeling. 

Av=-Av c [l — (rlr l ) 2 ], 0<r<rx 
and 

&v=Av m n<Cr<ir G . 

As an over-all balance is maintained, 

(Av)rdr= — (Av)rdr, 

and therefore eq 22, neglecting the D 2 /V a term, 
and assuming uniform lateral difTusivity, can be 
put in the form 



h=(2/D l V a rl) 



r \Av)rdr\ dr/r+ 

(Av)rdr dr/r 



This is integrated after substitution of the above 
expressions for Av, and the result is 

h=(2/L> 1 V a r%) 

{ (A^r*/4)[ln ^-3/4 + {njr,) 2 - 

(r 1 /r ) 4 /4]+(19/768)A^M). 

When the diameter of the channel is small com- 
pared with the tube diameter, such that r o ^>10r! 
and Av c r = 2Avmrl, this expression becomes, 
approximately, 

h = 0MAv m r,IV yiValD^^r.ln-O.l), (30) 

where the velocities and difTusivity are super- 
ficial and apply to the packing. 

With a central hole having an inside diameter 
from one-tenth to one-hundredth that of the tube, 
In r /r\ — 0.7 will lie between 1.6 and 3.9. On 
comparing this case with case IV, it is seen that 
for the same Av m in the packing, meaning that the 
total channel flow in milliliters per second is also 
the same, the h for axial core channeling is from 
8 to 20 times that for symmetrical wall channeling. 
Although this result may seem surprising, it is to 
be remembered that the axial channel presents 
far less area to the packing through which the 
compensating effect of cross diffusion can occur. 

From Poiseuille's law, the mean liquid velocity 
in the channel relative to the mean liquid velocity 
in the average pore within the packing is propor- 
tional to (r l jr v ) 2 , where r p is the equivalent pore 
radius. Hence, for a bead packing with 38 per- 
cent of void space, 



V a +Avj2 



V a l.38-Ax„ 



^(n/r,) 2 , 



or, approximately, 

Avj2V a =(AvJV a )(r /r i y=2,63(r l /r p y-l, 

when AVJV a <0.1. 

The equivalent pore radius for a bead packing 
has been determined experimentally to be 0.165X 
bead diameter [2]. 

Table 3 gives the values of h computed from eq 
30 for a 2-cm inside diameter tube packed with 
beads of 100-mesh per inch, r p = 0.0042 cm, 
assuming 1.0 KC\ at 18° C and a voltage gradient 
of 10 v/cm.__For these conditions, F a =0.00498 
cm/sec and Di = 1. 55 X10~ 5 cm 2 /sec. 
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Table 3. — Core channeling in 2-cm inside diameter tube, 
100-mesh beads, 10 v/cm 



Core channel radius 

Equivalent pore radius 

n/r* 


Core channel 
radius, r\ 


h 


cm 
4 


cm 
0. 0168 
.0252 
. 0336 
.0420 


cm 

0.068 

1.6 
15 
92 


6- 


8 


10 





Table 3 shows that when a single continuous 
core channel has reached a size about five times 
the mean pore size of the packing, the resultant 
h will be about 1.0 cm. In such a case isotope 
separation between ends of a 10-cm length of tube 
would be negligible. This result assumes, how- 
ever, that the radial diffusion coefficient has the 
value indicated, this being the minimum value 
corresponding to molecular diffusion. Actually, 
when the tube is set horizontally, M will have a 
much larger effective value than assumed because 
of thermal convections in planes normal to the 
axis. Hence, the values of h in table 3 should be 
regarded as maximum. They show merely that 
very poor results can be due to a single relatively 
small core channel. It is for this reason that bead 
packings should be held under pressure. 

Case VI. Face-Cooled Diaphragm, NonUniform 
Porosity. — The diaphragm of uniform porosity 
was considered under case II. When the dia- 
phragm is not uniformly porous, various regions 
can be found in which the porosity is above or 
below a mean porosity. Two factors of equal 
importance determine the performance of the 
packing, the first being the magnitude of the 
velocity departure due to the porosity departure; 
the second is the lateral distance across the face 
of the diaphragm between a velocity crest and a 
velocity trough. 

The general case does not permit a rigorous 
treatment. However, the order of magnitude of 
h can be gotten by assuming an equivalent dia- 
phragm in which the regions of porosity departure 
are parallel, the mean velocity departure and 
the mean crest-trough spacing of the actual case 
and the equivalent case being the same. 

The profile of the equivalent case is shown in 
figure 7 and is expressed analytically by 

kv=—kv m cos (irx/a). 



Equation 22a is applicable. As lateral diffusion 
vanishes at limits a and 0, as indicated in figure 5, 
the integration should be made between these 
limits. The result is 



A=0.0505(Az; m a/F a ) 2 y a /M+ D % \ '„, 



(31) 



this being the h for each symmetrical pair of 
strips, and therefore the h for the entire diaphragm. 
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Figure 7. — Profile of a diaphragm with variable porosity. 

From this result it is apparent that small values 
of h can be achieved only when the product, 
(velocity departure) times (crest-trough spacing), 
is small. 

Consider two diaphragms wherein the velocity 
departures are 10 percent of V a in one and 2 
percent of V a in the other. Assume a voltage 
gradient for each of 10 v/cm at 18° C, so that 
VJD l = 322 and 5 2 /V a == 0.0031 for 1.0 N K(1 
solution. The resulting values of h are given in 
table 4 for various crest-trough spacings, a. 

From the results in table 4 it can be seen that a 
nonideal diaphragm will enable appreciable isotope 
separation even though the centers of the regions 
of porosity departure are laterally a fraction of a 
centimeter apart, but, if several centimeters apart, 
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Table 4.— Diaphragms with variations in porosity, 10 v\cm 



dvJVa 


a 


h 


8V m /V a 


a 


h 




cm 


cm 




cm 


cm 


0.10 


0.5 


0.044 


0.02 


0.5 


0. 0047 


.10 


1.0 


.17 


.02 


1.0 


.0096 


.10 


2.0 


.65 


.02 


2.0 


.029 








.02 
.02 


4.0 
8.0 


.12 

.48 















only when the maximum velocity departure is less 
than about 2 percent of the countervelocity, V a . 
Case VII. Isothermal Capillary Film. — Consider 
the unobstructed passage formed between two 
elongated parallel walls that are closely spaced 
apart a distance 2w, are extremely smooth, and 
strictly parallel. If temperature effects be neg- 
lected, the velocity departure in terms of the 
thickness coordinate will be 

A^=(F a /2[3(z/w) 2 -l]. 
Equation 22a is applicable to this case and gives 



h=(2/105)V a w 2 ID+D/V a: 



(32) 



a result that is nearly identical to that for the 
capillary tube, case /. 

In practice, the walls should be good heat 
conductors for external cooling, and the inner 
surfaces, at least, should be electrically non- 
conductive. The film will run hot at the mid- 
plane, but the fact that viscosity times conductance 
is not precisely constant will affect the numerical 
coefficient in equation 32 by only a few percent 
and can be disregarded. Calculations for h can 
be based on values of T T a and D at 18° C. 

Thermal convection can be used to advantage 
to avoid a breadthwise velocity departure by 
arranging the walls so that length is horizontal 
and breadth vertical. A value of h calculated as 
in case 777 should be added to eq 32 when a 
breadthwise departure is present. The effective 
diffusion coefficient breadthwise can be found 
from eq 12 and 22a, 

In a second arrangement the film can be con- 
fined by vertical concentric cylinders. In this 
case, one of the cylinders should be rotated to 
avoid circumferential variations in composition 
and temperature. By proper selection of spacing 
the profile would be analogous to figure 2, the 
theory permitting lower values of h than indicated 
by eq 32. 



IV. Production Conditions 

Economical production of enriched material 
depends primarily on the values of e and h dis- 
cussed in sections II and III of RP1765 [2]. The 
larger the e and the smaller the h, the lesser will 
be the power requirement per unit rate of pro- 
duction. In general, the method is applicable, 
though not necessarily more economical than 
other methods, when e exceeds unity by as little 
as 0.1 percent and when h is less than 0.1 cm. 
In the following discussion a specific example will 
be worked through for e= 1.0030 and A=0.020 cm. 

The simplest procedure for continuous produc- 
tion is to maintain constant anode composition by 
exchange with an unlimited supply of normal K. 
After a build-up period, material enriched in 39 K 
can then be very slowly bled from the cathode. 
Another procedure, analogous to continuous frac- 
tional distillation, is to feed normal material to a 
point intermediate to the ends and to withdraw 
from both anode and cathode compartments. In 
this case the apparatus provides a stripping sec- 
tion yielding material at the anode enriched in 
41 K. The enriching section for 39 K may be treated 
as though an anode were located at the point of 
feed and is thus equivalent to the first procedure. 
The following treatment will be focused on the 
39 K enrichment. 

1. Build-up Period 

In bringing the cathode to the desired enrich- 
ment, the procedure is to operate at zero sampling 
rate (total reflux) while maintaining at the cathode 
end an exact balance between counterflow of solu- 
tion and K ion migration. Other factors remain- 
ing constant, it follows that the countervelocity 
must slightly increase as the enrichment increases; 
thus the maximum fractional increase in going 
from normal abundance to 100 percent 39 K would 
be (e— 1)/15.2, which is the difference between the 
mean ion velocity of normal K and the ion velocity 
of pure 39 K. This increase is well within the 
limits of the control and of no actual importance. 
However, assuming a uniform normality and a 
precise balance at the cathode, it follows that the 
countervelocity will be too high at the anode end 
by the same fraction. The result is a slight wash- 
ling out of salt, the voltage gradient increasing 
slightly toward the anode, and the tube thus 
balances itself throughout. More important is 
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the tendency of increased temperature to increase 
ion speed and destroy the balance, but here also 
the balance throughout the tube is restored by a 
change in normality of salt and voltage gradient. 
Hence it is strictly correct to suppose an exact 
mean balance for all sections of the tube, but it is 
not correct to assume uniform normality of salt 
except when temperature effects at the ends of 
the tube may be neglected. 

The velocity balance at any section requires 
v l c l + v 2 c 2 =V a {c l -\-c 2 ) or v 1 —V a =(v 1 — v 2 )c 2 /(c 1 +c 2 ). 
The net transport of light isotope across any sec- 
tionis (vi— V«)ci— D'dci/dsor (vi—v 2 )ciC 2 l(ci+c 2 ) — 
D'dci/dz mole cm" 2 sec -1 . Each factor in the 
foregoing is an average for the cross section, dis- 
tant 2 from the anode. 

The exact differential equation governing the 
total reflux build-up of concentration as a function 
of time and distance, and applying to a tube 
having nonuniform temperatures lengthwise, is 

bjD'dd/bz) d(v 1 -v 2 )c 1 c 2 /(c 1 +c 2 ) = bc , ) 

dz 2 bz ll ' K } 

In a very long packed tube it would be permis- 
sible to neglect temperature nonuniformities at 
the ends. In this case the effective dill'usivity, 
D'j the difference in the velocity of the ionic iso- 
topes, Vi — v 2 , and total concentration, Ci + c 2 , can 
be considered invariant. Equation 33 in terms 
of mole fractions then takes the form 

VNtlbsf-^bNiNdbz^YpbNi/ttt, (34) 

which involves the close approximation h= 
D'/V a ^D'lv 2 . 

Attempts to obtain an exact solution of eq 34 
in order to obtain the dependence of Ni on distance 
at the time the cathode has reached a desired 
enrichment have not been successful. This solu- 
tion is necessary to obtain the exact preliminary 
holdup of 39 K and the exact time at which pro- 
duction can be commenced. However, as it would 
apply only to a long tube, and since relatively short 
tubes or diaphragms are ample for 39 K concentra- 
tion, this failure is of secondary importance to the 
present problem. 

Returning to eq 33, which applies to a short 
tube, integration over the equivalent length 
V c /A of the cathode volume gives 

bcifbz=(vi—V2)c 1 c 2 /D'(ci+C2)- 

(VJAD^dcJdt^l, (35) 



showing that the concentration gradient at the 
front of the cathode volume is less than its 
maximum possible value at the concentrations 
by the amount (VJAD^dcJdt. This fact will be 
considered in the subsequent approximation for 
the build-up time of a short packed tube. 

2. Production Equilibrium 

When the cathode reaches the desired enrich- 
ment, withdrawal is begun at such a rate as to 
maintain the enrichment constant. Thereafter, 
the concentration distribution along the tube 
undergoes a further change until a steady dis- 
tribution is reached. 

Let AV a be the departure from the "total 
reflux" countervelocity, V a , this departure being 
due to continuous sampling at the cathode at a 
rate holding the cathode enrichment constant. 
When a steady distribution is reached a transport 
balance taken across the tube at distance z from 
the anode gives for the light component 

(vi--V a +AV a )c 1 --D'dc 1 /dz=AV a p 1 - eath , (36) 

and for the heavy component, 

(vi-V a +AV a )c 2 -D'dc 2 /dz=AV a c 2 _ cath , (36a) 

where the c's and D' refer to mean values for the 
cross section, and c x - caih) c 2 - calh refer to the moles 
per cubic centimeter of the two components in the 
product. 

Equation 36 is divided by c { and eq 36a by c 2> 
the latter is subtracted from the former, and the 
terms arranged to give 



dz~- 



D f {dcilci—dc 2 lc 2 ) CiC 2 



~CiC 2 (v 1 —V 2 )—AV a (C2C 1 - eath — CiC 2 -cath) 



(37) 



To proceed further, the approximation that 
temperature does not depend on z is assumed. 
The normality of salt is invariant and eq 37 can 
be written in terms of mole fractions as 



dz=- 



RhdN } 



-i?(6-l)A^ 1 2 +[i?(6-l) + lW 1 -iVl- 



(38) 



where v l -v 2 =v 2 (e-l)g*V a (e'-l),h=D'V a *!D'lv2, 
and R=V a /AV a g±v 2 /AV a . The integral of eq 38 
is had by a well-known formula. The result is 



^ftgJ ^aM + ft- Vg . 2aN 1 - mod 8 +b+J q 



hR V 2aN x 
Vg n_2aiVi 



+ J + V# ZaN^anoae + b— ^ 



} 



(39) 
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where a= -R(e-l), b = R(e-l) + l, c= -N^ cath , 
and <[=b 2 —4ac. Equation 39 gives the relation 
between distance from anode, z, and mole fraction 
light, Nij when the distribution has reached a 
steady state after prolonged production at reflux 
ratio, R. 

Example. — Consider the production of material 
having an abundance ratio greater than normal 
by 50 percent; ™K/ 41 K=N 1 /N 2 =21.3 at the cath- 
ode. Provided sufficient quantities can be pro- 
duced, this enrichment would suffice for many 
purposes such as experiments for tracing potassium 
metabolism. It is desired to know (1) the length 
of packing and reflux ratio for maximum efficiency, 
(2) the efficiency of the yield in grams of total 
enriched potassium per kilowatt-hour, (3) the 
required build-up time under total reflux operation. 

The case chosen is for a bead-packed tube, 1 
cm inside diameter with anode abundance main- 
tained normal at 14.2. The cathode volume, 
V c is 1.0 ml, and the equivalent length of cathode 
volume based on a packing of 38 percent voids is 
3.35 cm. The solution is taken as 1.0 N KC1. 
The temperature, normality, and voltage gradient 
will be assumed invariant with distance for 
purposes of this approximate calculation. The 
applied voltage is chosen to be 10 v/cm; the current 
is 0.50 amp and the input is 5.0 w/cm (a higher 
input would result in boiling). 

The value of e, based on the consideration of 
section II of RP1765 [2], together with experi- 
mental data reported in other papers of this series, 
is assumed to be 1.003. This is a conservative 
value. 

The value of h is taken to be 0.020 cm, a value 
which appears conservative for the optimum con- 
ditions of operation discussed under case III of 
section III. 

Equation 39, with the above values of e and h, 
is shown plotted in figure 8, for various values of 
R and corresponding lengths of packing, I. The 
Rl versus I curve is also shown. 

As seen in figure 8, the minimum packing 
length, I, required to give the assumed enrichment 
is 2.7 cm, under total reflux equilibrium, i. e., 
E= oo and t= *> . The required reflux ratio for a 
packing of infinite length is 114. 

The power requirement is proportional to 
length of packing whereas the production in grams 
of enriched material per unit time is proportional 
to 1/R. Hence, maximum efficiency is had when 



Rl is a minimum. The Rl versus / curve of figure 
6 shows a fairly broad minimum at about 4.6 cm 
for which i?=170; these are the values of packing 
length and reflux ratio which will give the most 
efficient yield although the curve shows that 
considerable latitude is permissible in the direction 
of greater packing length and lower reflux ratio. 
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Figure 8. — Relation between reflux ratio, R, and length of 
packing, I, for fixed cathode enrichment. 

The total potassium arriving at the cathode 
when counterflow is zero is I^/F moles/sec, 
where I + is the total current times the transport 
number; the latter in this case is close to 0.50. 
At reflux ratio R, the production is {I + jFR moles) 
/sec. For the present example the yield is 



. , A 0.50X0.50 
yield = ^t-^t^ttt^a = 1-5 



1 0~ 8 mole -of 2 1.3 K/sec,= 



96,500X170 

5.5X10" 5 mole of 21.3 K/hr. 

The power input to the 4.6 cm-packing is 
4.6X5=23 w. By allowing 7 w for electrode 
losses, the power, exclusive of that for the control 
equipment, is about 0.030 kw. Hence 

yield efficiency = 5.5 X10" 5 /0.030-l. 8 X10" 3 

mole of 21.3K/kwhr=7.2X10- 2 gof 21.3 K/kwhr. 
The corresponding energy requirement is 14 kwhr/g 
of 21.3 K. The energy charges against the produc- 
tion are thus seen to be quite nominal. 

The solid curve, 1, of figure 9 shows the abun- 
dance ratio vs distance from the anode, s, as 
related by equation 39. This is the production 
equilibrium. Curve 2 is a rough sketch of the 
solution of eq 33 at the time the cathode reaches 
the production enrichment; this curve has a slope 
at z=l that is less than curve 1, as was seen in eq 
35. As at the time of curve 2 the build-up is still 



182 



Journal of Research 



diffusing toward the anode and is less near the 
anode than for curve 1, it follows that curve 2 
must cross curve 1. 
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Figure 9. — Distribution of abundance ratio: 1, production 
equilibrium; {R = 170, t=&>); 2 t end of build-up timer; 
(R = oo f t=13 hr); 3, total reflux equilibrium; (R — oo } 

t= oo). 

The difference between areas A and B of 
figure 9 represents the hold-up of 39 K yet to be 
accumulated in the packing after the production 
has begun. This difference is negligible in the 
present case because of the relatively large hold-up 
at the cathode. In consequence, the approximate 
hold-up for both curves 1 and 2 is given closely by 
graphical integration under curve 1. This proce- 



dure gives: hold-up of packing plus cathode volume 
equals 2.2 X10" 5 mole of 39 K in excess of normal 
abundance 

As the concentration gradient at the anode re- 
mains zero for a large part of the build-up time, it 
is a good approximation to assume the hold-up 
of 39 K to have been transported at the maximum 
rate /+ {N X N 2 ) anode (e-l)/F moles 39 K/sec, from 
the anode into the packing. This rate is in 
accordance with eq 9 of section I of RP1765 [2] 
and gives for the present example, 0.5X0.5 
(14.2/15.2) (1/15.2) X 0.003 / 96,500 = 4.8 X 10" 10 
mole 39 K/sec. Hence total-reflux build-up time 

= 2.2X10- 5 /4.8X10- 10 

= 4.6X10 4 see. 

= 13 hr. 
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